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ABSTRACT 

We present results from a deep Chandra observation of MS1137. 5-1-66, a distant (z=0.783) and massive 
cluster of galaxies. Only a few similarly massive clusters are currently known at such high redshifts; 
accordingly, this observation provides much-needed information on the dynamical state of these rare 
systems. The cluster appears both regular and symmetric in the X-ray image. However, our analysis of 
the spectral and spatial X-ray data in conjunction with interferometric Sunyaev-Zel'dovich effect data and 
published deep optical imaging suggests the cluster has a fairly complex structure. The angular diameter 
distance we calculate from the Chandra and Sunyaev-Zel'dovich effect data assuming an isothermal, 
spherically symmetric cluster implies a low value for the Hubble constant for which we explore possible 
explanations. 

Subject headings: galaxies: clusters: general — galaxies: clusters: individual 

(MS1137. 5-1-6625) — cosmology: observations — cosmological parameters — distance 
scale — X-rays: galaxies: clusters 



1. INTRODUCTION 

Rich galaxy clusters are useful cosmological probes. The 
richest clusters are thought to be massive enough to com- 
prise a fair sample of the Universe, i.e., their mass composi- 
tion should reflect the universal mass composition (White 
et al. 1993; Evrard 1997) and so provide an efficient lab- 
oratory for measuring the baryonic to dark matter ratio. 
The mere existence of massive galaxy clusters at high red- 
shifts can also place powerful constraints on the physical 
and cosmological parameters of structure formation mod- 
els (e.g., Peebles et al. 1989; Bahcall & Cen 1992; Luppino 
& Gioia 1995; Oukbir & Blanchard 1997; BahcaU & Fan 
1998; Donahue et al. 1998; Eke et al. 1998; Haiman et al. 
2001; Holder et al. 2001). The greatest leverage is pro- 
vided by the most massive and distant clusters; in fact, to 
constrain JIa one must use clusters with redshifts greater 
than 0.5 (e.g., Oukbir & Blanchard 1997; Holder et al. 
2001). For these reasons, we are particularly interested in 
observations of massive galaxy clusters at high redshift. 

A set of observations of distant clusters can be used 
to constrain the expansion rate and curvature of the uni- 
verse. The most noted approach is through the angular 
diameter distance relation. The theoretical value of the 
angular diameter distance Da{z) is cosmology-dependent; 
Da{z) can be calculated from measurements of the X-ray 
temperature and surface brightness profile of a cluster's in- 
tracluster medium (ICM) analyzed in conjunction with a 
measurement of the cluster's Sunyaev-Zel'dovich effect (c/. 
Birkinshaw et al. 1991; Birkinshaw & Hughes 1994; Myers 
et al. 1997; Hughes & Birkinshaw 1998; Reese et al. 2000; 
Patel et al. 2000; Grainge et al. 2002), a spectral distortion 



of the cosmic microwave background (CMB) radiation by 
the hot, ionized ICM (Sunyaev & Zel'dovich 1970, 1972). 

Another approach to measuring cosmological constants 
with clusters uses the idea that the fraction of the total 
cluster mass contained in this atmosphere, the gas mass 
fraction fg, traces the universal baryonic mass fraction, 
under the fair sample assumption. The fair sample as- 
sumption then implies also that the gas mass fraction in 
massive clusters should not vary. 

With a well-chosen sample which includes clusters at 
high-redshifts, both the above determinations of the clus- 
ter gas mass fraction (c/. Sasaki 1996; Pen 1997; Rines 
et al. 1999) and cluster distances can be used to constrain 
the geometry of the universe, since at high redshifts, the 
calculated cluster distances and gas mass fractions depend 
significantly on the cosmological parameters ^Im and JIa- 
Recent measurements oi Da{z) via the X-ray/SZE method 
(Mason et al. 2001; Reese et al. 2002; Jones et al. 2003) 
and of the gas mass fraction (Grego et al. 2001) in samples 
of clusters imply values of the Hubble constant and the 
curvature of the universe in good agreement with the val- 
ues determined independently with other methods. In the 
work of Reese et al. (2002), the angular diameter distance 
to 18 clusters, distributed widely in redshift, is calculated. 
For ilA=0.7. they find the mean Hubble constant 

for the sample is 60^4 ~t\^ km s^^ Mpc^^. For seven low- 
redshift clusters. Mason et al. (2001) find a mean Hubble 
constant of 66li[^±15 km s~^ Mpc~^ for this same cosmol- 
ogy. Jones et al. (2003) find a mean value of 65l7±15 km 
s~^ Mpc~^ for five moderate-redshift clusters. (The un- 
certainties, statistical and systematic, are reported at 68% 
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confidence.) These Hubble constant measurements agree 
within the uncertainties with the value derived from the 
Hubble Space Telescope Ho Key Project (c/. Mould et al. 
2000), 72±3 ± 7 km s'^ Mpc^i, and Type la supernovae 
(c/. Riess ct al. 1998), who find the Hubble constant to be 
~65, in the cosmology considered here. 

When determining cosmological parameters from clus- 
ter distances and mass fractions, one generally assumes 
the dominant cluster physics is well-represented by a sim- 
ple picture: galaxies and an ionized intra-clustcr medium 
in hydrostatic equilibrium in a large dark-matter poten- 
tial, supported by thermal pressure. In general, this is 
descriptive of clusters in the local universe, as evidenced 
by the strong correlations between X-ray luminosity and 
temperature (Mushotzky & Scharf 1997; Allen & Fabian 
1998; Arnaud & Evrard 1999; Donahue et al. 1999) and 
cluster size and temperature (Mohr et al. 2000); and by 
the uniformity of the gas mass fraction in massive clusters 
(c/. David ct al. 1995; Mohr et al. 1999). One can attempt 
to ameliorate the effects of departures from this simple 
picture on the accuracy of the results by determining cos- 
mological parameters from samples of clusters, in which 
some of these systematic effects will average out, and tak- 
ing care to choose the sample in an unbiased way. Indeed, 
much of the observational cluster work cited above (e.g.. 
Mason et al. 2001; Grego et al. 2001; Reese et al. 2002) 
uses such an approach. In this work, we investigate the 
validity of the hydrostatic, isothermal, splicric;ally sym- 
metric assumption by studying observations of a massive 
(and therefore relatively bright) distant cluster. 

We analyze a deep observation of the galaxy clus- 
ter MSI 137.5-^6625 taken with the Chandra Observa- 
tory's Advanced CCD Imaging Spectrometer (ACIS). 
MS 1137.5-1-6625 was observed as part of L. VanSpey- 
broeck's Cycle 1 Guaranteed Time Observations (GTO). 
These GTO observations, together with observations in 
subsequent Chandra observation cycles, form a collection 
of X-ray observations of galaxy clusters designed for mea- 
suring cosmological parameters to high accuracy. Over 
1,250 ks of Chandra time were scheduled for this project 
in observation Cycles 1 and 2, and it included over 40 clus- 
ters by the close of Cycle 2. The project was designed to 
include a sufficient number of clusters to overcome sys- 
tematic errors in the calculated cosmological parameters 
arising from cluster ellipticity and orientation (c/. Sulka- 
nen 1999). The clusters in the sample are massive (kT,. > 
5 keV) and distributed widely in redshift. The clusters 
were chosen on the basis of having high X-ray luminos- 
ity and/or high X-ray temperature. Those with extremely 
bright radio point sources at their centers were avoided. 
The MS 1137.5-1-6625 observation is one of the longest ob- 
servations scheduled for this project. 

The Chandra Observatory carries instruments partic- 
ularly well-suited to studying the distant, massive clus- 
ters in the cosmology project, owing to a combination 
of unprecedented angular resolution in the X-ray band; 
low background rates, particularly in the front-illuminated 
ACIS-I devices; and the ability to detect X-ray photons up 
to high energies (kTe --10 keV). The XMM X-ray observ- 
ing facility has substantially more collecting area, permit- 
ting investigation of large-scale temperature gradients in 
distant clusters with observations of moderate length. The 



Chandra facility, however, has significantly better angular 
resolution than XMM, and as is discussed in this paper, 
this becomes quite important for distant clusters; their 
compactness does require the higher angular resolution of 
Chandra to accurately model the cluster's spectrum and 
spatial structure. The two facilities are therefore comple- 
mentary for study of high-redshift clusters. 

The Chandra observation reveals that while the X-ray 
image suggests this massive, distant cluster is relaxed, 
symmetric, and spherical, its story is actually consider- 
ably more complicated. We analyze the Chandra data 
in conjunction with SZE data and compare it with op- 
tical data. Wc determine it has a remarkably compact 
density distribution, find evidence of asymmetric temper- 
ature structure, and measure an angular diameter distance 
which implies a very low value for the Hubble constant. 

In this paper, we review previous observations of 
MS 1137.5+6625 in Section 2 and describe the Chandra ob- 
servations and the ICM temperature and density profiles 
we derive from them in Section 3. In Section 4, we derive 
the ICM cooling time, ICM mass and total mass profiles 
from the results of Section 3, and calculate the cluster dis- 
tance and infer the Hubble constant from the X-ray and 
SZE data in Section 5. We discuss these results, evaluate 
the cluster's suitability for cosmological work, and suggest 
possible interpretations in Section 6. 

2. PREVIOUS OBSERVATIONS 

MS1137. 5-1-6625 was serendipitously discovered in the 
Einstein Extended Medium Sensitivity Survey (EMSS) 
(Gioia et al. 1990) and was subsequently identified as a 
distant galaxy cluster (Henry et al. 1992; Gioia & Lup- 
pino 1994). It was confirmed by Donahue et al. (1999) to 
be at a redshift of z=0.784. Assuming f^M =0.3, riA=0.7, 
ft.=0.65 (which we will use throughout this paper), this red- 
shift implies an angular diameter distance (Da) of 1655.50 
Hq^ Mpc and a scale of 8.0 kpc arcsec"^. 

Further observations of the cluster confirmed that it is 
massive and compact. From a 70 ks ASCA observation, 
Donahue et al. (1999) find the cluster's best-fit rest-frame 
emission-weighted intracluster medium temperature to be 
kTe = 5.7ti;i keV, at 90% confidence. From Keck II spec- 
tra of 22 cluster members, they ascertain a velocity disper- 
sion of 884ti[^| km s'^ Clowe et al. (2000) obtained deep 
optical J?-band images of MS 1137.5-1-6625 with the Keck 
II 10 m telescope and /-band images with the University of 
Hawaii 2.2 m telescope. The optical light of the cluster has 
a compact distribution, and the centroid of this distribu- 
tion is coincident with the brightest cluster galaxy (BCG). 
There are two giant strong gravitational lensing arcs, at 
5.5" and 18.0" from the BCG, and several other arc candi- 
dates, indicating a high surface mass density. Clowe ct al. 
(2000) performed a weak lensing analysis of the data, and 
found the mass distribution is as compact as the light dis- 
tribution. Using aperture densitometry, they found the 
minimum cluster mass within a radius of 2' centered on 
the BCG to be (4.9 ± 1.6) x lO^^'/igs^Mg. A second mass 
peak to the north of the cluster may or may not be asso- 
ciated with the cluster, but only contributes about 10% of 
the total mass. 

3. Chandra observations 
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MS 1137.5+6625 was observed from September 30 to 
October 2, 1999 for an elapsed live time of 117.712 ks. The 
aimpoint of the ACIS-I3 chip was ~ 1' south and 2.5' east 
of the cluster center. We selected events with the stan- 
dard ASCA grades and rejected data taken during times 
with background rates more than 20% greater than the 
mean. The effective exposure time after this filtering was 
107.045 ks. In Figure 1, we present a raw image, without 
a background subtracted, in the 0.3 to 10.0 kcV energy 
range, blocked into 2" pixels. There are about 4000 clus- 
ter photons in the image. The background rate in the map 
is about 0.15 counts arcsec"^; the average emission rate 
over the cluster area is about four times brighter. In the 
image, the cluster looks relaxed and symmetric. 

3.1. Temperature Structure 

To determine the cluster's emission-weighted temper- 
ature, we extract the cluster spectrum within 63" (506 
^65 ^P'^) from the cluster's center. We estimate the back- 
ground in the same region using a number of long inte- 
grations on blank fields compiled by Markevitch. (See 
http://hea-www.harvard.edu/~maxim/axaf/acisbg/ for a 
complete description of the background datasets and their 
reduction.) We use the background data set taken co- 
evally with the MS 1137.5+6625 observations, so that the 
Chandra focal plane was at the same temperature (—110 
C). We use Vikhlinin's ca/crm/ package to derive response 
matrices and effective area files for the cluster spectrum; 
this software ensures that the variation of the energy re- 
sponse across the 13 chip is taken into account. Before 
extracting the spectra, we mask out the pointlike sources 
in the observed field. This has little effect on the num- 
ber of cluster photons we can retrieve; only one source, 
-^45" southwest of the cluster center, was within the clus- 
ter spectrum extraction region. The spectrum contains 
~4200 cluster photons and ~1800 background photons. 

We fit this spectrum, grouped into channels with a min- 
imum of 20 counts to ensure applicability of the fit 
statistic, to a single-temperature MEKAL plasma model 
in the 0.7-7.5 keV energy range, using XSPEC. The lower 
energy bound is chosen so that we keep only the best- 
calibrated data. We make our best estimate of the back- 
ground spectrum by an iterative process of varying the 
background spectrum normalization and noting how this 
affects the fit statistic for the best-fit model. Since the 
hard and soft components of the X-ray background orig- 
inate in different physical processes (cosmic X-rays dom- 
inate the background at energies below ~5 keV and cos- 
mic rays dominate above ~5 keV), the background in the 
MS 1137.5+6625 observation may depart from the blank- 
sky background differently at low and high energies. To 
mitigate these effects while keeping the high-energy cluster 
photons, we choose an upper bound of 7.5 keV for our fit. 
We find the best estimate of the background to be 4.5% 
higher than the blank-sky background file for our epoch. 

We find the best-fit emission-weighted mean tempera- 
ture for the intra-cluster medium to be 6.64+J;|q at 90% 

confidence. The best-fit metal abundance is 0.26lg'2i 
times the solar abundance as determined by Anders & 
Grevesse (1989), again at 90% confidence. The cluster 
spectrum, its best-fit model, and the fit residuals are 
shown in Figure 2 along with the two-parameter confidence 



regions for temperature and metal abundance. The neu- 
tral hydrogen column density is not strongly constrained 
by the data in this energy range. The 90% confidence in- 
terval is Uh = 2.2812^28 ^ 10^°cm~^, consistent with the 
Galactic value for this direction: 1.14x 10^°cm~^ (Dickey 
& Lockman 1990). The reduced statistic for this fit is 
1.051 for 150 degrees of freedom. 

Since the density profile suggests that the gas in the 
cluster center may be dense enough to cool (as discussed in 
Section 4), we also fit the cluster's spectrum to a MEKAL 
plasma plus a MEKAL cooling flow component. We fix 
the lower cutoff temperature to 0.2 keV and tie the upper 
temperature of the cooling flow to the MEKAL plasma 
temperature. The best fit cooling flow upper temperature 
is consistent with the emission-weighted temperature we 
found above, the cooling flow mass deposition rate is con- 
sistent with zero, and the reduced x^ is not improved. We 
conclude that there is no spectral evidence for a cooling 
flow component in the cluster. 

We compare these results to those obtained by Donahue 
et al. (1998, 1999), who found kTe = 5.7tll keV and an 

iron abundance of 0.431^:3? times solar abundance, at 90% 
confidence. These results are consistent with the Chandra 
results. Some difference may be expected, as there are 
a number of bright point sources in this field which can 
be distinguished by Chandra and their emission removed 
from the cluster spectrum, but which would likely con- 
taminate the ASCA cluster spectrum. The point sources 
in the 13 chip have a background-subtracted sum of about 
3800 counts in the 0.7-10.0 keV band. The brightest of 
these pointlike sources, located ~ 60" to the northeast of 
the cluster center, contributes ~ 1000 of the 3800 counts. 
This emphasizes the importance of good angular resolu- 
tion for studies of high-redshift clusters; it is difficult to 
determine accurately the ICM temperature and luminos- 
ity of high-redshift clusters without the ability to exclude 
point sources. 

To investigate the presence of structure in the tempera- 
ture or metallicity distribution in this cluster, we make a 
hardness ratio map using the adaptive binning technique 
described in Sanders & Fabian (2001). This technique en- 
sures a constant signal-to-noise ratio across the map; we 
set the fractional error in each bin to be 0.10. We estimate 
the background in each band from the MS 1137.5+6625 
dataset. In Figure 3 we show a map of the cluster in the 
2.0-5.0 keV band divided by the 0.43-2.0 keV band with 
the surface brightness contours of Figure 1 overlaid. This 
ratio is sensitive to temperature variations, moderately 
sensitive to metal abundance variations, and essentially in- 
sensitive to neutral hydrogen absorption variations. There 
appear to be significant fluctuations, with the center of the 
cluster softer than the outer regions. The hardness ratio 
image appears more complicated than the raw image, sug- 
gesting variations in temperature and/or abundance may 
not be present in the ICM density. Such behavior is not 
uncommon in clusters at low redshifts. Donnelly et al. 
(2002) derive temperature maps for a complete sample of 
nearby clusters observed with the ASCA facility. Many of 
these clusters which appear regular in surface brightness, 
even in images with much better photon statistics, exhibit 
quite complicated temperature structure, suggesting that 
temperature and abundance variations, likely created dur- 
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ing cluster merging, arc cither created frequently or persist 
for a long time. 

We can quantify the structure in this map further by 
identifying the hardness ratios with those expected from 
a MEKAL plasma. A MEKAL plasma at our best-fit 
emission- weighted temperature of 6.7 keV and metallic- 
ity of 0.25 solar gives a hardness ratio of -^0.28 in these 
bands. For a hardness ratio of 0.35, a ~ 7 keV plasma 
must have an abundance of ~ 1.0; or, if the abundance is 
to be 0.25, the plasma must be at a temperature of -^10 
keV. 

We attempt to confirm this temperature variation by 

fitting spectra extracted from different regions. The first 
is the region within 17" and the second is the annulus be- 
tween 17" and 47.5". The best fit temperature and abun- 
dance for the inner region are Q-9lLi\l ^^'^ ^-"^^-o.li 
times solar, at 90% confidence. The best fit temperature 
and abundance for the outer region arc 8.8^2 g ^cV and 
0.28lo!28 times solar. The data are consistent with a cooler 
inner region, but this is not statistically significant at 90% 
confidence. There arc not enough photons to support a 
temperature deprojcction analysis. 

3.2. Density Profile 

We parameterize the spatial distribution of the electron 
number density spherically-symmetric isothermal 

beta-model(c/. Cavaliere & Fusco-Femiano 1976, 1978). 
The beta-model has the form: 

/ ^2 \ -3/3/2 

ne(r) = rieo (^1 + , (1) 

where Ueo is the electron number density at the cluster's 
center; /3, though derived with a physical meaning, is in 
practice a fitted parameter; Vc is a characteristic length 
scale for the cluster, also a fitted parameter; and r is the 
radius from the center. Thermal bremsstrahlung X-ray 
emission from such a plasma results in a surface bright- 
ness profile of the form 

Sx{e) = Sxo[l + -^J , (2) 

where Sxa is the surface brightness at the center of the 
cluster, 9c is the cluster's projected core radius 9c = 
Tc/Da, and 9 is the angular distance from the cluster cen- 
ter. Using the CXC Sherpa package, we fit an image made 
in the energy range 0.3 keV to 5.0 keV to a generalized 
beta-model, in which two perpendicular axes can have dif- 
ferent core radii. (As noted by Fabricant et al. (1984), a 
biaxial ellipsoid with a beta-model electron-density distri- 
bution will project to a simple elliptical beta-model of this 
type for any orientation of the axes with respect to the 
observer.) The axis ratio, i.e., the ratio of the minor to 
the major axis, is 0.871q q7, at 90% confidence. We use the 
image analysis to quantify the deviation from circular sym- 
metry and to locate the cluster's center. However, using 
a biaxial or triaxial cluster model introduces a consider- 
able amount of ambiguity, cf. Grego et al. (2000). We use 
a spherically symmetric model, checking that the best-fit 
shape parameters from the two-dimensional fit are con- 
sistent with those from the spherically symmetric model, 
and leave the discussion of the effects of deviations from 
spherical symmetry until Section 6. 



Using the centroid from the two-dimensional fit (RA = 
11'*40"'22".4, Dec.= +66° 08' 14".7), we derive a one- 
dimensional surface brightness profile and fit it to the 
beta- model of Equation 2. We model the background (as 
a constant) concurrently with the cluster model, so that 
we can use the log-likelihood statistic (called the "Cash" 
statistic in Sherpa). This statistic, S = —2ln{L) where L 
is the Poisson likelihood, tends asymptotically to be dis- 
tributed like the statistic, so confidence intervals on 
the calculated quantity can be constructed by the devi- 
ation from the minimum statistic. Over the small angu- 
lar region within which we are fitting (within a radius of 
85") and in this energy band, the effective area of Chan- 
dra changes by less than 3% and the background (as esti- 
mated from the Markevitch background files) is constant. 
We also allow for an offset between the centroid we deter- 
mined from the image analysis and the surface brightness 
peak. In Figure 4 we show the data, the best fit model, 
and the fit residuals; the best-fit model has /3=0.63 and 
a core radius of 13.3", or 106 h^^ kpc. The confidence 
intervals on the surface-brightness shape parameters are 
shown in the right panel of the figure. For the confidence 
interval calculation, the offset between the centroid and 
peak is fixed to its best-fit value of 2.1", reasonable for a 
cluster with the observed axis ratio. If we require the peak 
to coincide with the centroid, the best-fit core-radius de- 
creases to 10.9", the best-fit (} is 0.56, and the fit statistic 
degrades by AS'=1.6. The small core radius of this clus- 
ter reinforces the importance of Chandra's high angular 
resoution for studying high-redshift clusters. 

We fit this profile again, after eliminating the central 
region, to see if the small (and marginally statistically sig- 
nificant) central excess is affecting the shape parameter fit, 
and we find no appreciable change in the fitted parameters. 

We recover the central gas density from the normaliza- 
tion of the fitted XSPEC model to the spectrum. The 
XSPEC normalization is as follows: 

^^"= 4.(i+.)^i?i y 

where the integral is taken over the volume projected onto 
the area within which the spectrum was extracted, here, 

within a circular region with radius R. All variables are 
in cgs units. The central density can be found from this 
relation: 




where C, is the distance from the cluster center along the 
line of sight in units of the core radius. 

To determine the uncertainty in the central electron den- 
sity measurement due to the spatial fits, we calculate 
at a large number of parameter pairs (grid points) dis- 
tributed over the 90% confidence interval on P of Fig- 
ure 4b. Determining the confidence interval on n^o from 
the relevant confidence area for 9c and /3 is necessary to 
correctly account for the strong correlation between pa- 
rameters. We add in quadrature the spectrum normal- 
ization uncertainty and the uncertainty from the spatial 
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fits. We find a central density of 0.0136^ 



hO.0016 -3 
-0.0013 "-"^ 



In 



Figure 5, we illustrate the relative contributions of the un- 
certainties; we plot with a solid line the ne{r) calculated 
for each grid point within the 90% confidence interval (the 
contribution from the spatial fit) and with a dashed line 
the 90% confidence interval on ne(r) which includes the 
contribution from the uncertainty in the XSPEC normal- 
ization. We note that uncertainty is dominated by the 
uncertainty in the shape parameters. 

4. CALCULATED QUANTITIES 

4.1. Cooling Profile 

We calculate the cooling time of the ICM as a function 
of cluster radius, using the ratio (5fcTe/2)/(e//(Te)ne(r)), 
where e//(Te) is the free- free emissivity of a plasma with 
our best-fit abundance (line emission contributes very lit- 
tle cooling at 7 kcV) . The cooling time in Gyr as a function 
of cluster radius is shown in the left panel of Figure 6. In 
the right panel, we compare the cooling time to the age 
of the universe at the cluster's epoch, which is 7.53 hg^ 
Gyr in our chosen cosmology. The cooling time of the 
ICM appears to be shorter than the cluster's ago at radii 
below 25-70 h^^ kpc. However, we do not sec convinc- 
ing evidence for a cooling flow in the cluster's image or 
spectrum. For the cooling time in the center of the clus- 
ter to be longer than the cluster's age, one or more of our 
assumptions must be incorrect. There could be a source 
of heat to the gas at the cluster center (although there is 
no evidence for a radio-loud galaxy there). The cosmol- 
ogy we use could be incorrect, but for the cooling time to 
be longer than the cluster's age with certainty, h must be 
greater than 1.0 for an f2M=0.3 flat universe, or f2M=0.8 
if /i is ^ 0.65 in a flat universe. It is also possible that 
our simplified physical model for the cluster is incorrect; 
if the gravitational potential is steeper along the line of 
sight than it appears in projection, the central density we 
calculate would be overestimated, and the cooling time 
underestimated. A relatively modest change of OcQine of 
sight) /^^c (projected) = 0.8 can lower the central density 
by ^50%. Or the cluster may have formed or undergone 
a merger with enough power to disrupt cooling processes 
at some time more recent than tcool- 

4.2. Gas and Total Mass Profiles 

The gas mass profile follows directly from the electron 
number density profile (c/. Sarazin 1988): 



-ft 

Mg{R) = An J ne{r) Heinpr^ dr. 



(5) 



We use the best-fit abundance to calculate the mean 
atomic mass firUp and assume /U is constant throughout 
the gas, and that the plasma is fully ionized. We calculate 
the cluster's total mass under the assumption that the 
ICM is isothermal and in hydrostatic equilibrium in the 
cluster potential, with thermal support only. Under such 
assumptions, the total mass of a cluster within radius R is 

where G is the gravitational constant. 



The resulting gas mass fraction profile is shown in Fig- 
ure 5, where again we use solid lines to trace fg{r) for the 
allowed shape parameters and the dashed lines to delin- 
eate the 90% confidence interval which includes the un- 
certainties in Nxs a-iid Tg. We compare the left and right 
panels of Figure 5 and note that the statistical uncertainty 
is dominated by the uncertainty in Tg. 

We compare the gas mass fraction we calculate with the 
Chandra data with the gas mass fraction determined us- 
ing Sunyaciv-Zel'dovich effect data. Wc use the method 
described in Grego et al. (2001) which was used to de- 
termine the gas mass firaction from interferometric SZE 
data and published ASCA temperatures. The method 
uses the SZE data to constrain the cluster's shape pa- 
rameters P and 6c and X-ray spectra to determine the 
emission- weighted temperature, which is combined with 
the SZE normalization and shape parameters to deter- 
mine Ueo- To better compare the masses, wc redcrive 
the SZE fg using the Chandra temperature from Section 
3.1. The SZE/ C/ianrfro-determined gas mass within 63" is 
3.01°'q X lO^'^ft-^tf Mq and the gas mass fraction within 

63" is O.OGlQ Qg/igg^. The precision of the temperature 
measurement is greatly improved from that used in Grego 
et al. (2001); this is masked somewhat by our choice here 
to use the 90% confidence interval rather than the 68% 
confidence interval. We compare these to the gas mass 
and gas mass fraction determined in this work from the 
Chandra data alone: a gas mass 

Mq and a gas mass fraction of 0.157to;o44^65^^^- '^^^ 
suits are summarized in Table 1. The gas mass and gas 
mass fraction measurements are marginally consistent at 
90% confidence. 

To compare the total mass estimate with the weak lens- 
ing measurement, we calculate the total cluster mass pro- 
jected within a radius of 2'. We integrate the mass density 
along the line of sight to ~20 core radii, or about 1.5 Mpc; 
if we integrate out twice as far, the total mass changes 
by only ~15%. The projected cluster mass within 2' de- 
termined from the Chandra data is 6.0'^\'^xlO^'^MQhQ^ . 
This is quite consistent with the minimum projected mass 
of (4.9±1.6) X IO^^H^^Mq calculated from the Clowe et al. 
(2000) weak lensing estimate. 

We note here that the cluster's lensing-derived mass pro- 
file appears more compact than expected. In their weak 
lensing analysis, Clowe et al. (2000) compare the calcu- 
lated mass profile to the Navarro et al. (1996) universal 
profile and find MS 1137.5-1-6625 is more centrally concen- 
trated than is usual for massive clusters in this scenario. 
(The Navarro et al. (1996) profile has since been noticed 
to generally predict steeper density distributions than ob- 
served, so this is particularly striking.) MS 1137.5-1-6625 
also exhibits a number of lensed arcs in the cluster cen- 
ter, again indicating a high central mass concentration. 
Clowe et al. (2000) find these points suggestive, but not 
conclusive, of an elongation of the cluster along the line 
of sight. Our analysis of this cluster shows a core radius 
of ~ 100 kpc. This is certainly more compact than 
the bulk of nearby clusters, whose core radii are generally 
around 150 h^^ kpc (Jones & Forman 1999; Mohr et al. 
1995). Only 11 of the 45 clusters in the Mohr et al. (1995) 
sample have core radii as small or smaller than this, yet 
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the temperature of MS 1137.5+6625 has a higher emission- 
weighted temperature, and is therefore presumably more 
massive, than 28 of the clusters in the sample. The clus- 
ter may then fall into the compact-core radius group of 
Ota & Mitsuda (2002), which divides clusters into two dis- 
tinct groups, one with core radii around 60 kpc, and one 
with core radii 220 kpc. Although a full investigation into 
correlations between core radius size and other properties 
has yet to be done, Ota & Mitsuda (2002) find evidence 
that clustc^rs with c;D galaxies have small core radii and 
that the luminosity-temperature relationships for the two 
cluster populations have different normalizations, though 
similar slopes. 

5. CLUSTER DISTANCE 

We derive the cluster's angular diameter distance by 
comparing the cluster's X-ray emission and Sunyaev- 
Zel'dovich Effect (SZE) data. The method closely follows 
the method used in Reese et al. (2000), except that in this 
work we use Chandra imaging and spectral data instead 
of i?05^T imaging and ASCA spectral data. 

The SZE has a magnitude proportional to the Compton 
y-parameter, i.e., the total number of scatterers, weighted 
by their associated temperature, 

I ne{l)Te{l)dl, (7) 

where k is Boltzmann's constant, ctt is the Thomson scat- 
tering cross section, nie is the electron mass. For an 
isothermal beta-model, the projected SZE effect would be: 

AT = f^^^^^^TcMBDA /rfC^rne^ = ATo (l + J)^ 

' (8) 

where AT is the thermodynamic SZE temperature decre- 
ment/increment and f(x,Te) is the frequency dependence of 
the SZE, with x = hv/kTcMB, where Tcmb is the tem- 
perature of the cosmic microwave radiation Tcmb=2.728 
(Fixsen et al. 1996). In the non-relativistic and Rayleigh- 
Jeans limits, ,f(x.T,) ~ —2; we apply the relativistic cor- 
rections of Itoh et al. (1998) to fifth order in kTf,/mf,(?, so 
that = -1.907. 

Extending Equation 2, the cluster's projected X-ray sur- 
face brightness is 

= 4.(1 + .)^ "^^^^"-"^^-^ = + 

(9) 

where Sx is the X-ray surface brightness in erg s ^ cm ^ 
arcmin"^, z is the cluster's redshift, j\h is the hydrogen 
number density of the ICM, = K^niTe, abimdancc) is 
the X-ray cooling function of the ICM in the cluster rest 
frame in erg cm"^ s^^ integrated over the redshifted Chan- 
dra band, and 5x0 is the X-ray surface brightness at the 
center of the cluster. We convert the X-ray detector counts 
to cgs units using the Chandra response matrices avail- 
able from http://asc.harvard.edu/cal/Links/Acis/acis/; 
the conversion factor is 7.91x10"^^ erg cm"^ . One 
can solve for the angular diameter distance by eliminating 

We perform a joint fit to the interferometric SZE data 
and the Chandra image, in which the X-ray pointlike 
sources have been masked out. We compare the X-ray 



data to a beta model with a constant cosmic background; 
this model is exposure-corrected before the fit. The loga- 
rithm of the Poisson likelihood is then calculated. The in- 
terferometric SZE observations provide constraints in the 
Fourier {u-v) plane, so we perform our model fitting in the 
u-v plane, where the noise properties of the data and the 
spatial filtering of the interferometer are well defined. The 
Gaussian likelihood is calculated using the SZE model and 
the SZE data. There are no detected point sources in the 
SZE data. We use a downhill simplex method (Press et al. 
1992) to search the parameter space and maximize the 
joint likelihood of the cluster position, /?, 6^, Sxo, ATo, 
and cosmic X-ray background. 

Each data set is independent, and likelihoods from each 
data set can simply be multiplied together to construct 
the joint likelihood. From this we can generate confidence 
regions for the parameters. 

Uncertainties in the angular diameter distance from the 
fit parameters are calculated by gridding (stepping) in the 
interesting parameters to explore the AS likelihood space. 
As a compromise between precision and computation time, 
we grid in ATq, Sxo, and 9c allowing the X-ray back- 
ground to float while fixing the positions of the cluster 
(both SZE and X-ray). From this four dimensional AS 
hyper-surface, we construct the 90% confidence interval 
for Da due to Sxo, ATq, /3, and 6c jointly. The corre- 
lations between the beta-model parameters require this 
treatment to determine accurately the uncertainty in Da 
from the fitted parameters. We emphasize that these un- 
certainties are meaningful only within the context of the 
i/3)/2spherical isothermal f3 model. 

, In the Reese et al. (2002) work, the calculated angular 
diameter distance of MS 1137.5+6625 using ROSAT and 
ASCA data was SUdtHilh'^ Mpc (statistical uncertainty 
only), at 68% confidence. In our analysis here, with a more 
accurate emission- weighted temperature and improved X- 
ray imaging data, we find an angular diameter distance 
(at 90% confidence) of 3439t^^^^/i-^ Mpc, which implies 
a Hubble constant of 3lli[y km Mpc"^ (with statistical 
uncertainties only) . The Chandra observation moderately 
improves the precision of the measurement, but do not 
rectify it with the Reese et al. (2002) sample; the derived 
Hubble constant is still ~ 2.5cr from the 18 cluster sam- 
ple average.. MS 1137.5+6625 is a peculiar cluster, and in 
the next section we discuss our calculations and attempt 
to discern how this cluster may deviate from a relaxed, 
spherical, isothermal cluster. 

6. DISCUSSION 

The X-ray and optical images of the distant galaxy clus- 
ter MS1137. 5+6625 show the X-ray surface brightness and 
optical light have a regular, symmetric, and quite compact 
distribution. However, several lines of evidence suggest 
that the cluster structure is more complicated: 

• The hardness ratio map and spectral analysis sug- 
gest that the cluster has a relatively more compli- 
cated temperature and metal abundance structure 
than density structure. 

• The gas mass fraction calculated from the 
SZE/ Chandra data is only marginally consistent 
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with the gas mass fraction calculated from the X- 
ray data. 

• The Hubble constant calculated for this cluster is 

low by a factor of ^2. (This is closely related to the 
discrepancy in the gas mass fraction and gas mass.) 

These calculated quantities arc summarized in Table 1. 
The most compelling argument that MS 1137.5+6625 devi- 
ates from isothermal, spherical symmetry is that the Hub- 
ble constant value derived from it is so low. The X-ray 
and SZE data sets are not sensitive enough individually 
to allow us to definitively test our assumptions, but taken 
together they do provide some insight. Here, we discuss a 
number of possible explanations for our observations. 

6.1. SZE Data Systematics 

The possible contribution to systematic imcertainty to 
the Hubble constant and gas mass fraction measurement 
by the SZE data itself is very low. It is possible that there 
are undetected radio-bright point sources in the cluster 
field. Unless they are unfortunately placed on the side- 
lobes, point sources will decrease the measured SZE effect 
and increase the measured Ho- Combining the SZE ob- 
servations with 1.4 GHz NVSS (Condon et al. 1998) ob- 
servations suggests that possible undetected point sources 
in the MS 1137.5+6625 field have a negligible effect on the 
derived Hubble parameter (Reese ct al. 2002). We find 
no evidence for a radio halo in this cluster in the NVSS 
image. Radio haloes, i.e., large scale diffuse radio emis- 
sion, would again depress the measured SZE and therefore 
would cause an overestimate of the Hubble constant. 

Anisotropics of the CMB of other types which might af- 
fect the SZE data are at very low levels. A limit to the 
primary anisotropics at these angular scales was set by 
Dawson ct al. (2001) and Holzapfel et al. (2000). They 
find the Rayleigh- Jeans temperature of these fiuctuations 
to be less than 14 /xK, and this in turn should contribute 
less than 2% uncertainty to Ho- Confusion from the kine- 
matic SZE effect (the spectral distortion of the CMB made 
by the cluster moving with respect to the CMB rest frame) 
should also be quite small; for an 8 keV cluster moving 
with a 300 km s~^ velocity component in the line of sight 
direction, the kinetic SZE effect will be ^4% of the ther- 
mal SZE (and the sign of the effect will depend on whether 
the cluster is moving towards or away from the observer) 
leading to an 8% error in the Hubble constant. It is there- 
fore quite unlikely for the SZE to be severely in error due 
to peculiar velocity of the cluster gas. 

6.2. Geometry 

The Hubble constant is derived under the assumption 
that the cluster can be characterized by a spherically sym- 
metric distribution. This will be in error by the ratio of 
the projected cluster size Lpnoj to the line of sight cluster 
size Llos'- Hmeas = Htrue X Lproj/Llos- We have mea- 
sured the apparent axis ratio to be 0.87tQ;oy for the ellip- 
tical beta model; the core radius measured by the circular 
beta model is roughly the geometric mean of the two axes. 
For an apparently circular distribution, the line of sight 
axis could be longer or shorter than the projected axis, 
leading to a prolate or oblate three dimensional cluster, 



respectively. For an apparently elliptical mass distribu- 
tion, the true three-dimensional distribution has another 
degree of complexity in addition to the degree of prolate- 
ness or oblateness; the axis of symmetry could have any 
angle of inclination with respect to the plane of the sky, 
allowing the symmetry axis an essentially infinite range 
of lengths which would produce the same projected im- 
age. Additionally, there is no a priori reason that clusters 
need to be ellipsoidally symmetric; in fact, for axis ratios 
of less than l/\/2, the concentrically ellipsoidal isothermal 
beta- model for the density distribution is unphysical (c/. 
Grego et al. 2000) if the the cluster is indeed in hydrostatic 
equilibrium. The gravitational potential supporting such a 
beta-model requires negative dark matter density. A more 
sophisticated model for the non-spherical density distri- 
bution is clearly required in such cases, but is beyond the 
scope of this paper. We note that for measuring cosmo- 
logical parameters from an ensemble of clusters, increasing 
the complexity of the model is not strictly necessary. For 
example, (Sulkanen 1999) shows that from a sample of tri- 
axial clusters, one can reconstruct an unbiased estimator 
for the Hubble constant by using a spherical model. 

In the simplest picture, the error in the Hubble constant 
is approximately the ratio of the apparent cluster size (i.e., 
the geometric mean of the two core radii) to the line of 
sight size. In such a picture, for MS 1137.5+6625's mea- 
sured Hubble constant to equal the sample mean of Reese 
et al. (2002), L (projected) /L (line of sight) nmst be 0.58. 
As a fiducial point, in a fiux limited sample of clusters 
fitted with elliptical beta-models, Mohr et al. (1995) mea- 
sure an axis ratio for the sample of 0.8;*]q ]^^. For the low 
Hubble constant to be due solely to an elongated cluster, 
it would have to be have two unusual characteristics for 
a cluster chosen from an unbiased sample: it would have 
a much longer axis ratio than generally observed in pro- 
jection, and would have its symmetry axis nearly aligned 
with the line of sight, basically a long cigar pointed at the 
observer. If the symmetry axis were not nearly along the 
line of sight, the true axis ratio would need to be even 
greater. 

Since such an elongation along the line of sight is 
so unexpected, this may indicate that the criteria by 
which MS 1137.5+6625 was detected favored such clusters. 
MS 1137.5+6625 was originally identified in the EMSS 
(Gioia et al. 1990), using a 2. '4 x 2. '4 detection cell, a 
size which would just fit the cluster emission we observe 
with Chandra. It was detected with a signal-to-noise ratio 
of 5.3 in a survey with minimum S/N for detection of 4.0. 
MS 1137.5+6625 is the coolest and least luminous of the z 
> 0.5 clusters in the EMSS survey. If there does exist a 
population of similarly luminous, highly elongated clusters 
at high redshifts, it is likely that only those aligned very 
nearly along the line of sight would be detected by this 
survey; a cluster also near the detection cutoff but which 
was elongated in the plane of the sky would deposit less 
of its flux into a single detection cell and may not be de- 
tected. For example, if MS 1137.5+6625 had an intrinsic 
axis ratio of ^0.75 but with the major axis in the plane 
of the sky, it is unlikely that it would have been detected 
with sufficient significance to be included in the survey. 

The selection eff'ects from the original point-source al- 
gorithm used for the EMSS were discussed in Pesce et al. 
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(1990), and reprocessing of the Einstein data led Lewis 
et al. (2002) to eonclude that the EMSS survey systemat- 
ically missed clusters of low surface brightness. Similarly, 
simulations of the EMSS detection algorithm by Ebeling 
et al. (2000) suggest the algorithm may miss clusters with 
significant substructure. 

6.3. Temperature Gradients 

Although we have no compelling evidence for a cool- 
ing flow at the cluster's center, we do have evidence of 
structure in the temperature and abundance distribution. 
Temperature structure could affect both the derived shape 
parameters and the derived emission-weighted tempera- 
ture. 

Temperature structure will affect the X-ray and SZE 
data differently. It should not strongly affect the X-ray- 
derived spatial shape parameters; for a given emission 
measure, the emitted flux in the band we use for spatial 
fitting only changes by ^ 20% for gas with kTg from 4 keV 
to 10 keV. Also, for the X-ray spatial fitting, the data are 
azimuthally averaged, so that such structure would have 
little effect. 

The SZE measurements, by virtue of the spatial filtering 
inherent in interferometry, are insensitive to temperature 
gradients on scales larger than a few arcminutes, but can 
be quite sensitive to smaller scale gradients. As the SZE 
effect is proportional to x Tg, the density inferred at 
any particular point will be over (under) estimated from 
an under (over) estimate of the temperature. Tempera- 
ture structure could result in the beta-model fitted to the 
SZE differing from that fitted to the X-ray data. We do 
not expect this poor fit to the beta-model to strongly af- 
fect the calculated Hubble constant, because the joint-fit 
to the model is dominated by the X-ray imaging data in 
this case. It may, however, affect strongly the gas mass 
fraction determined from the SZE data, as this is derived 
from the SZE spatial fits. 

If MS 1137.54-6625 had a non-hydrostatic pressure 
structure, then we may expect to see evidence of this in 
the SZE data. The low surface brightness of the SZE-effect 
makes it difficult to image on small scales, however, and 
with the current data, we are unable to make a conclu- 
sion about non-hydrostatic pressure. Future generations 
of SZE instruments may have such capability. 

The temperature structure of the cluster could con- 
tribute to the low value of the Hubble constant in an- 
other way. The temperature we use for the Ho calcula- 
tion is the emission-weighted temperature, and the rele- 
vant quantity for the calculation is actually the density- 
weighted temperature. As a fiducial point for the discus- 
sion, we compare the emission-weighted temperature and 
the density-weighted temperature for a cluster with the 
following structure: /3 equal to the best-fit value from our 
model fit, and gas within two core radii (~200 kpc) at 
half the temperature of the gas from two core radii to ten: 
Te(r < 2rc) = 0.5 x Te(2rc < r < lOrc). Temperature 
structure of this magnitude, while unusually pronounced, 
has been observed by Chandra in the cluster Abell 1835 
(Schmidt et al. 2001), a massive cooling-fiow cluster with 
highly luminous optical emission-line nebulae. In such a 
scheme, the emission-weighted temperature would be 0.80 
times the density-weighted temperature, and the Hubble 



constant would be 0.64 times the true Hubble constant. 
For the underestimation of the Hubble constant to be 
due solely to an underestimate of the temperature, the 
emission-weighted temperature must be ~0. 65-0. 70 times 
the density-weighted temperature. We stress that this is 
a simple calculation; the exact effect depends also on the 
specifics of the instruments used to make the measure- 
ments. However, even if MS 1137.5-1-6625 had as extreme 
a temperature gradient as the atypical Abell 1835, this 
could not explain the entire Hubble constant underesti- 
mate. 

7. SUMMARY 

In conclusion, the massive, distant galaxy cluster 
MS 1137.5-1-6625 appears compact and regular in its op- 
tical and X-ray images, but evidence suggests that it is 
not an isothermal, spherical cluster. The Hubble constant 
calculated for the cluster from X-ray imaging and spectra 
and interferometric SZE data is unusually low. 

Many of the known systematic effects in the X-ray/SZE 
Hubble constant measurement tend to make the Hubble 
constant higher, and few effects make it lower. (A thor- 
ough discussion of these systematic effects can be found 
in Reese et al. (2002).) Two effects which can make the 
measured Hubble constant artificially low are elongation 
along the line of sight and temperature structure in the 
ICM. We have evidence suggestive of a line of sight elon- 
gation. The optical lensing data indicate the surface mass 
density is quite high, as evidenced by strongly lensed arcs, 
and that the mass distribution is very compact. The X- 
ray data show the ICM density distribution is unusually 
compact, as well. We see structure in the hardness ratio 
map which could be interpreted as temperature structure, 
with cooler or less metal-enriched gas at the center of the 
cluster. 

We also compare other calculated quantities besides the 
Hubble constant to investigate the systematic effects. The 
SZE-derived gas mass is completely insensitive to geome- 
try, though it is quite sensitive to temperature structure; 
if the emission-weighted temperature were biased low, the 
SZE gas mass would be an overestimate of the true gas 
mass. The X-ray-derived gas mass is essentially insensitive 
to temperature structure, but has a moderate dependence 
on geometry; if the cluster were elongated along the line of 
sight, the X-ray gas mass would be a slight overestimate. 
The X-ray-derived gas mass is 2.ll°;} h^^^"^ x lO^^^o and 
the SZE-derived mass is S.Oti^ h^^ x 10"A/o. Within the 
uncertainty the masses agree and we cannot place strong 
limits on temperature structure and elongation with the 
current datasets. 

We conclude with two points: one, that cluster appear- 
ance is a poor criterion for choosing sources for cosmolog- 
ical tests. MS 1137.5-1-6625 appeared as a possibly ideal 
cluster for cosmology, with its regular, compact, round, X- 
ray and optical images, and no evidence of a strong cooling 
flow. However, it yields a Hubble constant measurement 
which is strikingly low. We stress that a cluster sample 
must be selected via well-defined, objective criteria. And 
two, that selection effects in the parent sample from which 
the clusters are drawn should be well- understood, i.e., al- 
though for cosmological work we generally select clusters 
from existing X-ray surveys on the basis of luminosity 



rather than surface brightness, any selection for surface 
brightness in the X-ray survey itself will still be present in 
the luminosity-selected sample. This situation can be ame- 
liorated by future cluster surveys, such as the SZE survey 
discussed in Holder & Carlstrom (2001), which have selec- 
tion functions different from the X-ray and optical cluster 
surveys. 
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Table 1 

Summary of Derived Properties 



Data 



gas mass total mass 

within 63" projected within 120" 



within 63" 



(for n 



M 



Da 
(Mpc) 

= 0.3, = 0.7) 



Chandra 
SZE 

SZE/Chandra 
Keck Weak Lensing 



-5/2 
65 

+ 1.0 , -2 
-1.0 "-65 



2 h 
■^■^-0.1 "65 

3.0 



6.o: 



-1.7 
-1.2 



4 q+1-6 



0.157 



+0.022 1,-3/2 
-0.044 "-65 



0.06 



-0.07 . -1 
-0.03 "-65 



0.00+1854,-1 
O'±o»_2087"'65 




RA (J2000) 

Fig. 1. — Raw 0.3 keV to 10.0 keV image of MS11H7. -3+6625, blocked into 2" pixels. There are approximately 4,000 cluster counts in 107 ks 
of time. No background subtraction has been performed. The contours are logarithmically spaced from 0.55 to 30 counts, where the map has 
been smoothed over 4 pixels. The aimpoint is ~ 1 arcmin south and ~2.5 arcmin east of the cluster's center; the instrument PSF is nearly 
unchanged from the aimpoint at the cluster's position. 
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Fig. 2. — Left panel: the Chandra ACIS-I spectrum, best-fit model, and residuals for the full cluster spectrum for MS 1137.5+6625. Right 
panel: the two-parameter confidence region for emission-weighted temperature (kTe) and abundance. The first, third, and fourth regions 
correspond to 68.3%, 90%, and 99.0% confidence on the two parameters jointly. The projection of the second region onto either axis marks 
the 90% confidence interval on the single parameter. The marks the best-fit value. 
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Fig. 3. — Hardness-ratio map of MS1137. 5+6625. The cluster image in the 2.5-5.0 keV energy range is divided by the 0.5-2.5 keV image. 
The background is subtracted from each image and the fractional error is a constant value of 0.15 across the map. An emission- weighted 
temperature and abundance of 7.0 keV and 0.245 solar would give a hardness ratio of ~ 0.165. Lower values correspond to lower temperatures 
and/or metallicity. The surface-brightness contours from Figure 1 are overlaid. 
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Fig. 4. — Left panel: One-dimensional surface brightness profile of cluster MS1137.5+6625 with background, its best-fit beta-model, and fit 
residuals in units of the standard deviation <t. . Right panel: Two-parameter confidence region for j3 and Oc for the one-dimensional beta- model 
fit to MS1137.5+6625. The confidence contWe use values Qm = 0.3 and 0.7.ours correspond to 68.3%, 90%, 93% and 99.0% confidence 
on the two parameters. The projection of the 93% confidence interval onto an axis gives the single-parameter 90% confidence interval. The 
best-fit model is marked with a '+'. 
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Fig. 5. Left. Electron number density versus radius in kpc. Right. The gas mass fraction of the cluster as a function of cluster radius in 
kpc. The solid lines show ne(r) and fg(r) for parameter grid points within the 90% confidence region of the fitted iriodcl. The dashed lines 
show the full 90% confidence intervals for ne(r) and fg(r), which includes the uncertainty in the spectrum's normalization for ne(r), and in 
the spectrum's normalization and emission-weighted Te for fg{r). The calculations are done using values CIm = 0.3 and CIa= 0.7. 
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Fig. 6. — Left panel: The cooling time of the ICM in Gyr as a function of cluster radius in kpc. Right panel: The cooling time of the ICM 
in units of the age of the universe at the cluster's redshift, as a function of cluster radius in kpc. The calculations are done using values 
= 0.3 and Qa= 0.7. Again, the solid lines show the cooling times calculated for the parameter grid points within the 90% confidence region 
of the fitted model and the dashed lines show the full 90% confidence intervals including the uncertainty in Nxs and Tg. The age of the 
universe at the cluster's epoch is 7.53 Gyr hg^ . 



